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Abstract 

Using the Korringa-Kohn-Rostoker coherent-potential approximation in the atomic-sphere 
approximation (KKR-ASA CPA) we have studied the changes in the electronic structure 
and the superconducting transition temperature Tc in Nbi^xB2 alloys as a function of x. 
We find that the variation in the electronic structure of Nbi-xB2 alloys as a function of 
X is consistent with the rigid-band model. However, the variation of T^, obtained using 
the Allen-Dynes equation within the Gaspari-Gyorffy formalism to estimate the electron- 
phonon matrix elements, does not follow the expected trend. We associate this disagree- 
ment to the use of a constant Wrms in the Allen-Dynes equation over the whole range of 
vacancy concentration, thereby indicating the importance of lattice dynamical effects in 
these systems. 
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1 INTRODUCTION 



The discovery of superconductivity in M gB2 and the fact that the carbides and the 
nitrides of Nb show superconductivity have renewed interest in finding supercon- 
ductivity in the borides of Nb, in particular, the diborides. However, the search for 
superconductivity in NbB2 and A^6i?2-based alloys is not new. Earlier work by Zei- 
gler et al [1] and Hulm et al [2] showed B deficient NbB2, namely NbBi g^ to be 
superconducting with superconducting transition temperature, Tc, of around 1.2 K, 
while Schriber et al [3] found NbB^ (synthesized under pressure) with x ~ 2 to be 
superconducting with Tc of 9.2 K. 

Recently, Yamamoto et al [4] and Kotegawa et al [5] have reported superconduc- 
tivity in Nbi-xB2 alloys when synthesized under pressure for x ranging from to 
0.48. They find that Nbi_xB2 samples become superconducting for x = 0.04, and 
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the Tc increases to around 9.2 for x = 0.24, thereafter Tc remains essentially 
unchanged up to x = 0.48. 

In an effort to theoretically understand the changes in the electronic structure of 
Nhi-xB^ and NbBi g4 alloys and their superconducting properties as vacancies are 
introduced either in the Nb plane or the B plane, we have carried out first-principles 
electronic structure calculations of these alloys as a function of vacancy concen- 
tration. We have used Korringa-Kohn-Rostoker coherent-potential approximation 
[6,7] in the atomic-sphere approximation (KKR-ASA CPA) method for taking into 
account the effects of disorder, Gaspari-Gyorffy formalism [8] for calculating the 
electron-phonon coupling constant A, and Allen-Dynes equation [9] for calculating 
Tc in Nbi-^B2 and NbBi,g4^ alloys. 

Here, we like to point out that the recent theoretical work of Shein et al [10] to 
model disordered NbQ,'j^B2 alloy using a 12-atom supercell may not be adequate. 
This can be seen as follows. Each metal atom in AlB2-tygQ diborides has D^h 
symmetry and thus has \2 B neighbors at the vertices of a hexagonal prism. In 
addition, the central metal atom also sees 8 others through the faces of the B12 
prism at a distance comparable with the metallic diameter. Thus a minimum of 
24-atom supercell would have been more appropriate. However, their conclusion 
that the superconductivity in A^fei_^i?2 alloys is governed more by the changes that 
take place in the lattice properties when vacancies are introduced in the metal plane 
of the diborides, is consistent with our results. Before describing our results, we 
provide some of the computational details of the present approach. 



2 Computational Details 

2.1 The densities of states 

The charge self-consistent electronic structure of A^6i_xi?2 and NhBi ^s alloys has 
been calculated using the KKR-ASA CPA method with x ranging from to 0.48. 
The lattice constants for Nbi_xB2 alloys are taken from the experimental work 
of Yamamoto et a/[4]. To determine the electronic structure of NbB2 using the 
present approach, we have used the lattice constants a and c from our earlier full- 
potential linear muffin-tin orbital calculations[l 1]. The absence of Nb and B atoms 
in the Nb and B planes, respectively, were modelled by introducing empty spheres 
at those sites. The disorder was accounted for by the coherent-potential approx- 
imation [12,13,14], which has been very successful in describing reliably many 
physical properties of disordered alloys [7]. The exchange-correlation potential 
was parametrized as suggested by Perdew-Wang within the generalized gradient 
approximation [15,16]. The Brillouin zone (BZ) integration was carried out with 
1215 k-points in the irreducible part of the BZ. For densities of states (DOS) calcu- 
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lations, we added a small imaginary component of 1 mRy to the energy and used 
4900 k-points in the irreducible part of the BZ. For valence states the angular mo- 
mentum cut-off was kept as l^ax = 3, and the core states were recalculated. For 
Nh and B the ratio of the atomic sphere radius to the Wigner-Seitz radius was kept 
fixed at 1.294 and 0.747, respectively. The empty spheres inherited the radii of the 
corresponding atomic spheres. The present calculations were carried out using the 
scalar-relativistic Schroedinger equation. The Green's function is calculated in a 
complex plane with a 20 point Gaussian quadrature. 



2.2 The electron-phonon coupling and the superconducting transition tempera- 
ture 



The superconducting transition temperature Tc was calculated with the Allen-Dynes 
equation. 



Ic = — exp 
" 1.2 ^ 



1.04(1 + A) 



A-yu*(l + 0.62A) 



(1) 



In our calculation, the value of the phonon frequency uin for NbB2 was taken 
from Ref. [11], fi* was set equal to 0.09, and the electron-phonon coupling con- 
stant A was calculated using the Gaspari-Gyorffy formalism with the charge self- 
consistent potentials of Nbi-xB2 and NbBi^g-^ alloys obtained with the KKR-ASA 
CPA method. This entails writing A in terms of Hopfield parameter, rj, as 



where the sum is over the basis atoms in the primitive cell, M is the atomic mass 
and < u"^ >^/^ is an average phonon frequency. The average phonon frequency 
< o;^ >^/^ is taken to be the root-mean- square frequency ujrms —< co'^ >^^'^. In the 
Gaspari-Gyorffy formalism the spherically averaged Hopfield parameter is given 
by 

''"^^(2/ + l)(2/ + 3) N{Ep) • 



The total DOS at the Fermi energy E^, N{Ef), and the /-resolved DOS, Ni{Ef) 
as well as the electron-phonon matrix elements, 

Mi,i+, = -<j>i(EF)(t>i+i(EF) \(Di{Ef) - 1)(Di+,(Ef) + 1 + 2) + (Ef- V{S))S^] (4) 
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are calculated from the charge self-consistent potentials obtained with the KKR- 
ASA CPA method. In Eq. (4), <j)i{Ep) is the amplitude of the /— th partial wave at 
the sphere boundary, S, evaluated at Ep, Di{Ep) is the corresponding logarithmic 
derivative, and V{S) is the one-electron potential at -S". 



3 Results and Discussion 



In the following we describe the results of our calculations for NhB2 and its alloys 
in terms of (i) the densities of states, (ii) the electron-phonon coupling and (iii) 
the superconducting transition temperature. We first describe how the densities of 
states, the electron-phonon coupling and the superconducting transition tempera- 
ture change for low concentrations of vacancies in either the Nh or the B plane. 
Then, we describe the changes in these properties as a function of increasing va- 
cancy concentration. 



3.1 The densities of states of NhB2, Nbos9B2 and NbBig^ 



3.1.1 NbB2 

In Fig. 1 we show the calculated total DOS and the partial DOS together with the 
/-resolved partial DOS of NbB2. The overall structure of the DOS is consistent 
with the previous work[10,17]. The calculations based on the tight-binding LMTO 
method yields a value of total DOS at the Fermi energy equal to 4.65 states/Ry- 
cell [17], while the corresponding value from the full-potential LMTO calculations 
is found to be 4.28 states/Ry-cell. The present calculations, which is based on the 
KKR-ASA method, gives N{Ef) to be equal to 4.41 states/Ry-cell. At Ep the 
contributions of Nb and B to the total DOS is found to be 3.44 and 0.97 states/Ry- 
cell, respectively. 



The characteristic features in the DOS of NbB2, as shown in Fig. 1, are {i) a low 
lying broad peak around 0.90 — 0.75 Ry below Ep due to B s— and Nbd— states, 
(ii) another broad peak at around 0.45 — 0.25 Ry comprising of Bp—, Nbp— 
and Nbd— states, (Hi) a relatively constant DOS from 0.7 — 0.5 Ry primarily of 
Bp— and Nbd— states, and (iv) a pseudo-gap at 0.14 Ry, which appears in many 
transition-metal diborides [17]. Thus the dominant presence of Nb d— states at Ep 
in NbB2 makes it quite different from MgB2 which has largely B p— states at Ep. 
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Fig. 1. (a) The total (solid line) and the partial densities of states of Nb (dotted 
line) and B in NbB2, calculated with the KKR-ASA CPA method as described in 
the text, (b) The various symmetry-resolved (i-partial DOS of A^6 in NbB2 (dxy : 
dashed line, dy^ : dotted line and d3z2-i'. full line). Note that d^-i-yi and d^z are de- 
generate with dxy and dyz, respectively, (c) The symmetry-resolved p-partial DOS 
of B in NbB2 ipz '■ full line and '■ dotted line). In this case, Py is degenerate with 

Px- 

3.1.2 Nbo.99B2 

To understand how low concentrations of vacancies in the Nb plane of NbB2 would 
affect its electronic properties, we have used KKR-ASA CPA method to calculate 
the electronic structure of dilute Nbi-xB2 with ,r— 0.99 alloy. In this calculation 
the lattice constants were kept the same as that of NbB2. In Fig. 2, we show the 
calculated /-resolved partial DOS at the empty sphere site in A^6o.99^2- It is seen 
that upon 1% creation of vacancy in the Nb plane, the total DOS at Ep for Nbo,ggB2 
is reduced by ~4% relative to that of NbB2. We find that the creation of vacancies 
in the Nb plane do not lead to appreciable change in the B partial DOS at Ep. Thus 
the metal deficiency introduced in NbB2 enhances the ratio of B states to that of 
the total DOS at Ep. 




3.1.3 A^65i.93 

Motivated by the work of Otani et al on NbBi_gs alloy, which showed that a deple- 
tion of B in NbB2 led to an increase in the lattice parameter c from 6.16 to 6.25 
a.u. and a decrease in the lattice parameter a from 5.89 to 5.84 a.u., we have cal- 
culated the electronic structure of NbBi^g^ alloy. The calculation is carried out by 
disordering the B sub-lattice with empty spheres and keeping the Nb sub-lattice 
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Fig. 2. Top panel: The /-resolved partial densities of states (s : dot-dashed line, p 
: dashed line, d : solid line, and / : dotted line) of the empty sphere in A^6o.99-S2 
alloy, calculated with the KKR-ASA CPA method as described in the text. Lower 
panel: The symmetry-resolved d-partial DOS of the empty sphere in (dxy 
: full line, dy^ : dotted line and dsz'2_i: dashed line). Note that d^^^y^ and d^z are 
degenerate with d^y and dyz, respectively. 



ordered. The calculated densities of states at Ep are summarized in Table I, where 
we have also listed the lattice constants used in the present work. 

Our analysis of the DOS of NbBi,Qs shows that the creation of vacancies in the B 
plane largely reduces the d^^'^^i component (~12.5%) of the Nb rf-states and the 
B 6-states (~8.9%) at Ep. The creation of vacancies in the B sub lattice of NbB2 
thus leads to a major reduction in those states which are directed along the c— axis 
of the unit cell. 



Table I. The lattice constants a and c, in atomic units, used for NbB2 and NbBi_g3 alloys. 
The calculated partial (Nb and B) and total densities of states at Ep, in states/Ry-cell, are 
also given. 





a 


c 


N^\Ef) 


N^{Ep) 


N^°\Ep) 


NbB2 


5.847 


6.259 


0.497 


3.497 


4.491 


NbBim 


5.887 


6.161 


0.473 


3.294 


4.241 
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3.2 The electron-phonon coupling and the superconducting transition tempera- 
ture of NhB2, Nbo_ggB2 and NbBi g^ alloys 

The superconducting transition temperatures of NbB2, NhQ,ggB2 and NbBi,Q3 were 
calculated with the Allen-Dynes equation. The factor representing the Coulomb 
interaction fi* was 0.09. For Urms and coin, we used the values 464 cm^^ and 494 K, 
respectively. Our results for the Hopfield parameter r], electron-phonon coupling 
constant A and the superconducting transition temperature are summarized in Table 
n. From Table n, we see that r]B, the Hopfield parameter due to B, is relatively 
small in comparison to rfm, which reflects the fact that B contribution to the total 
density of states at Ep is small. Since we use the same ujrms for NbB2, A^6o.99-B2 
and NbBi^gs the changes in A are similar to that of r]. The calculated Tg for NbB2, 
Nbo.ggB2 and NbBi,gs are found to be 0.31K, 0.13K and 0.004K, respectively. In 
the case of NbBi_gs alloy, we find that the deficiency of B largely affects rjNb rather 
than TjB, which reduces Tc significantly. 

Table II. The calculated Hopfield parameters 77, in mRy/a\, the electron-phonon coupling 
constant A and the superconducting transition temperature Tc{K) for NbB2 and NbBi_g^ 
alloys. Note that Atot=AAr;,-i- 2\b for NbB2 and \tot=XNb+ l-93As for NbBi,g^. 





VNb 


VB 


^Nb 


Xb 


Atot 


Tc 


NbB2 


81.7 


I8.I 


0.053 


0.121 


0.294 


0.313 


NbBi,gs 


61.2 


13.9 


0.039 


0.089 


0.212 


0.004 



3.3 The densities of states of Nbi_xB2 alloys 

In Fig. 3, we show the DOS of Nbo.geB2, Nbo,84B2, A^6o.76-S2 and A^6o.52-B2 alloys 
using the KKR-ASA CPA method. The lattice constants for these alloys were taken 
from Ref. [4]. The changes in the DOS that arise in Nbi-xB2 alloys as a function of 
vacancies are (i) the inward movement of Ep with increasing vacancy concentra- 
tion, leading to a decrease in the total DOS at Ep for x < 0.16 and an increase for 
X > 0.16, and (^^) the broadening of peaks with increasing vacancy concentration 
due to disorder. We find that the movement of Ep follows the trend as suggested 
by the rigid-band model, and the DOS of NboMB2 alloy is similar to ZrB2. 

3.4 The electron-phonon coupling and the superconducting transition tempera- 
ture of Nbi_xB2 alloys 

From the resulting electronic structure we have calculated the variation in the Hop- 
field parameter, electron-phonon coupling constant and the superconducting tran- 
sition temperature of Nb-i-^B^ alloys as a function of vacancy concentration. Our 
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Fig. 3. The total (solid line) and the partial densities of states of Nb (dotted line) 
and B (dashed line) in (a) iV6o.96-B2, (b) (c) iV6o.76-B2 and (d) Nbo.52B2 

alloys, calculated with the KKR-ASA CPA method as described in the text. The 
vertical line represents the Fermi energy. 




results for 77 and A are summarized in Table ni. We find that both 77 and A de- 
crease with increasing vacancy concentration up to x ~ 0.25, thereafter they start 

to increase. In Fig. 4, we show the calculated as a function of x in Nbi^^B2 
alloys, together with the experimentally reported values. The parameters used in 
the calculation of Tc were jj,* = 0.09, Urms = 400 cm~^, and ujin = 494 The 
calculated trend in the variation of Tc is in disagreement with the experimental re- 
sults as shown in Fig. 4. In our opinion, the disagreement is very likely due to the 
use of a constant Urms value in the Allen-Dynes equation for the whole range of x 
from to 0.48 for Nbi-^B2 alloys. 

In the present approach, the calculation of Tc incorporates two important approxi- 
mations, (i) use of Gaspari-Gyorffy formalism for estimating the electron-phonon 
matrix elements, which is known to underestimate the strength of the coupling, 
and (ii) use of a constant phonon frequency over the whole range of x. Using 
Gaspari-Gyorffy formalism with concentration-weighted averages of the phonon 
frequencies, we have been able to describe reliably the trend in the variation of 
Tc in many MgB2-hased alloys [18,19,20]. Since the electron-phonon coupling in 
NbB2 is dominated by Nb and not by B, as is the case in MgB2, the Gaspari- 
Gyorffy formalism may be more suitable for NbB2 than for MgB2. This can be 
seen by comparing the T^. as calculated in the present approach with that of the 
linear response calculation of Ref. [11]. The lack of agreement between the present 
calculations and the observed Tc in Nbi-xB2 alloys as a function of x indicates the 
importance of incorporating the lattice-dynamical effects in these calculations. 
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Table III. The calculated variation of Hopfield parameters t], in mRy / a\ and the electron- 
phonon coupling constant A as a function of x in Nbi-xB2 alloys. Note that Xtot=i^- 
x)A7Vfe+2As for Nbi-xB2 alloys. 





TtATU 

'tlvb 


'IB 




Ar 


^tot 


NbBo 


80.3 


21.1 


0.070 


0.158 


0.387 




73.1 


18.2 


0.063 


0.136 


0.334 


jV Dfi no _Do 


65.2 


15.4 


0.057 


0.115 


0.283 


7V6o.84^2 


55.9 


12.4 


0.048 


0.092 


0.227 


Nbo.7eB2 


47.6 


12.9 


0.041 


0.097 


0.228 


^^0.68^2 


62.9 


21.0 


0.054 


0.157 


0.356 


NboQoB2 


73.3 


27.3 


0.064 


0.205 


0.454 


Nbo.52B2 


82.2 


32.1 


0.072 


0.241 


0.526 




Fig. 4. The calculated (filled triangle) and the observed (filled circle) superconduct- 
ing transition temperature Tc of 7V6i_^B2 alloys as a function of x. 

4 Conclusions 



We have studied the electronic structure of NbB2, NbBisi^ and Nbi^xB2 alloys 
using the KKR-ASA CPA method. We find that the vacancies on the metal plane 
reduce the number of Nh d-states at Ep. The changes in the density of states in 
Nh-i-xB2 alloys as a function of x reflect a rigid-band picture. The variation in Tc in 
Nhi^xB2 alloys as a function of x, calculated using the Gaspari-Gyorffy formalism 
disagrees with the observed trend. In our opinion, this disagreement arises due to 
the use of a constant Urms in the Allen-Dynes equation over the whole concentration 
range in our calculations. 
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